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Abstract 

The present paper develops the ideas introduced in cond-mat/ 03 12673. 
The construction of a hybrid discrete-continuous model of layered su- 
perconductors is briefly presented. The model bases on the classic 
Lawrence-Doniach scenario with admitting, however, long-range in- 
teractions between atomic planes. Moreover, apart from Josephson 
couplings they involve the proximity effects. The range of interac- 
tions, K, can, in principle, be arbitrary large. The solutions corre- 
sponding to the range K=2 are exposed. The fundamental excitations 
are understood as deviations from stable ground states. The formulae 
for energy of those excitations are constructed. The possible shapes 
of dispersion curves are analysed. For each type of shape the corre- 
sponding values of physically measurable quantities like effective mass 
and bandwidth are expressed by coupling parameters. 



1 Introduction 

The layered structure of superconductor implies a very strong structural 
anisotropy which may be characterised by relation between ^ the coherence 
length in the direction perpendicular to layers - and the interlayer distance 
s. In the limiting case s the 3D anisotropic Ginzburg-Landau model 
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is well applicable [1, 2, 3, 8]. On the other side, in the case < s, a much 
better description is given by Lawrence-Doniach model: a stack of 2D GL 
planes with Josephson links between adjacent planes [4, 9]. Similar model 
has been used in [5, 7, 6, 10]. However, there exist superconducting materials 
like fiig^rg CaCng (9g or Tl2Ba2CaCu2 0s, which have ^c/^ between 1.5 and 
2, i.e. too small for GL and to high for LD. 

The higher grade hybrid model (HM) [11, 12, 13, 14], presented here, 
patches the gap between those extreme cases, admitting Josephson-type cou- 
plings (we shall call them J-links) between distant planes. The range K of 
J-links is a global parameter of the model. K=0 denotes the stack of non- 
interacting GL planes, K=l - J-links between nearest planes only, K=2 - 
J-links between nearest and next nearest planes, and so on. 

Like in LD model we consider two continuous independent variables under- 
stood as in-plane coordinates x, y. The third independent variable is discrete: 
it denotes the ordinal number of the plane n. 

We shall use the following notation: 
ipn — the order parameter associated to the layer indexed by n 
ipn — its complex conjugate (c.c), 
niab - in-plane 



nic out-of-plane 



effective mass of superconducting current carriers. 



2 Free energy functional 

Contribution to the superconducting component of the free energy contains 
2 parts: from GL planes and from J-links. 

„ J Zniab z 

+|EE/ dxdy{CM' + \^nU')-lSn^n+,e-'P''- + c.c.)}, (1) 

neP qePn 

where P denotes the set of indices of all planes, Pn refers to planes which 
are J-linked to the plane n. Every J-link is represented by exactly one term. 
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The symbol D denotes the 2-dimensional operator of covariant derivative 

ic* 

Dp^dp- — Ap, p = x,y. (2) 
The third component of the vector potential apears in the exponent 

g* f{n+q)s 



Pqn 



p-^ r[n+q s 



3 Comparison with the anisotropic GL model 

To compare our hybrid model (HM) with the continuum GL model, we shall 
consider the infinite medium. In that case the component of free energy 
connected with J-links has the form 

Pj = ^ EE{2(C, - IqMnl' + 7.|V^n+,e-^- - ■ (4) 

n q 

For very weak fields Az and very slow dependence of ipn on n we have the 
correspondence rules which, in the long wave limit, allow us to pass from 
the hybrid to the 3D continuum. Applying the rules formulated in [14] one 
obtains 

^.-/ ^M^iDV'r+^?v7,p3V'r+aoiV'r+E[(c.-7.)iV'r+^/?iV'ii}. 

(5) 

It implies the modification of GL parameter ao to the form 



a 



= "0 + E(C</ - Iq) (6) 



and gives the following relation between the GL effective mass in 2;— direction 
and the coupling parameters. 

(^E?H) \ (7) 
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4 Field equations 

By computing the variation of the functional Tg with respect to i/Sn, one 
obtains the equations 

1 

-7, D^V^^ + aV'n + /?| V'nrV'n - 77 E Tgf^gnV'n+^e '^"^ = (8) 

where ^ 

a = "0 + 2 II ^9nCg (9) 

depends of n for finite P. The quantities 

« = {-A-,.... -2, -1.1, 2,.... if}. -,= {jli:;i3''^^ (10) 

have been introduced to describe the finite range of J-hnks. 

The expression for Joscphson current Ji describing tunelhng across the 
interplanar gap indexed by / (half integer if planes are indexed by integers) 
will have the form 

se* - 

= ^ E E {Tq'^qnV'nV'n+^e-'P''" - C.C.}, (11) 
qeQ nePig 

where 

Piq^{neP -.TKl <n + q}. (12) 



5 The ground states 

Consider now the plane-uniform states of HM in the absence of magnetic 
field. Let us concentrate on the case K — 2. The condition of vanishing 
Josephson current is equivalent to 

7l(^n^n+l - C.C.) + 72(^n^n+2 + Ipn-llpn+l " C.C.) = 0, (13) 

and the equations (8) take the form 

aipn + I3\lpn\^1pn - ^[7l(^n+l + 1pn-l) + 72(^n+2 + lpn-2)] = 0. (14) 
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Solutions with constant amplitude and difference of phase between adjacent 
atomic planes, obtained by the ansatz ipn = C e"*^ form 3 classes of states 
with = -a*/P: 

uniform : il^n = C, = + Ci + C2 - 7i - 72- 

alternating : ipn = (-1)"'C, a* = cto + Ci + C2 + 7i - 72, 

phase modulated : ipn — Ce^'^'^^ , ^ = arccos(— - — ). (15) 

472 

The phase modulared solutions exist if the parameters 71 and 72 fulfill the 
relation |7i| <4|72|. Then a* is connected with the coupling constants by 
the formula ^ 

= ao + Ci + C2 + 72(1 + A). (16) 

^72 

The straight lines 71 + 472 = 0, 71 — 472 = 0, and 71 = divide the 
plane (71,72) into 3 regions (shown in Fig. 1 of [14]). Both the uniform 
and the alternating solutions always exist. However, in the region : 
7i > 0, 7i + 472 > 0, only the uniform solution is stable, while in the 
region @ : 71 < 0, 71 — 472 < 0, only the alternating solution is stable. 
The region (S) : 72 < 0, 472 < 71 < —472, excludes the stability of both 
the uniform and the alternating solutions but, in contrast to that, ensures 
the existence and stabihty of the phase modulated solutions . In the sector 
(N): 72 > 0, I71I < 472, the phase modulated solutions exist but are unstable. 

Note that, irrespectively of the values and signs of 71, 72, a stable ground 
state solution always exists. 

For suitable relations between the coupling constants 71 and 72, one can 
make the parameter a* more negative than aQ. In consequence, the 3D 
superconductivity can turn out enhanced with respect to the 2D supercon- 
ductivity in the atomic planes. The presence of coupling constants (q allows 
to take into account the proximity effect between atomic planes. According 
to the formula (16) one can obtain a negative a* starting from a positive ckq. 
It means that under appropriate coupling one can obtain supercoductivity 
by stacking planes which originally are in normal state. 
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The HM with K=l has in general two couphng parameters. ( and 7. In 
particular cases it is related to known models: 
C = 7 > LD model, uniform solution, no enhancement. 
C = 7 < Theodorakis model, alternating solution, proximity effect, 
enhancement. 

6 The excited states 

The physical interpretation of the model requires adequate understanding of 
the excited states. In the present section we shall consider the excited states 
as deviations from the ground states determined in the previous section. 
The results presented here for K—2 may be directly generalized to the case 
of arbitrary K [13]. 

For two reasons the problem will be treated in the linear context: 

(a) The nonlinearity is present in the equations through the term I'^p only, 
hence the excitations which do not change the modulus [■^l (we shall call 
them phase excitations) may be treated in the framework of linear equations 
exactly. 

(b) For the excitations changing the modulus of the order parameter, one can 
solve the problem in linear approximation, valid for weakly excited states. 
From the microscopic point of view one speaks here on the collective excita- 
tions of the Cooper pair condensate. 

Let us first consider the phase excitations. We substitute the function 

i/jn — /e'^", where / = const. (17) 

into the equation 

(a + P\ij\'^)iJn - ^[jlii^n+l + i^n-l)i^n + 72(^n+l + ^n-l)V'n] = e^/'n- (18) 

Then we obtain 

(19) 

which should be independent of n. In consequence. On — n{9o + 9) + const 
and 

e{e) =a + f3f- [71 cos (^0 + ^) + 72 cos2(^o + 0)], (20) 
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where 6q corresponds to the ground state determined in the Section 5. The 
particular values 6*0 = and 9o = n represent the uniform and alternating 
ground states, respectively. The remaining values of represent the phase 
modulated ground states. 

The expression for the energy of a phase excitation has the form: 
e(^) - e(0) = 71 (cos ^0 - cos (^o + ^) +72(cos2^o - cos2(^o + 0)), (21) 
where e(0) denotes the energy of ground state 

e{0)^a{9o)+Pf. (22) 

Now we shall examine the modular excitations. Let us consider once more 
the equation (18) but this time we shall take into account that ■0n f+^i^n, 
where / denotes a uniform ground state. After linearization one obtains 
equation of the same shape, with in place of lipl"^ and Sip in place of ip. 
Taking excitation Sipn = ae-'^^^^ with arbitrary constant amplitude a and wave 
vector k (it is sufficient to consider excitations with kx = = ky) , after some 
calculations like in the case of phase modulated excitations, one obtains the 
excitation energy corresponding to wave vector k (with respect to the ground 
state A; = 0) 

€i — €k — €o — 7i(l — cos ks) + 72(1 — cos 2ks) (23) 

valid in the region of stability of uniform ground state: 71 > 0, 71 + 472 > 0. 
Similar procedure perfomed for the alternating ground state delivers the for- 
mula 

£2 = -71(1 - COS ks) + 72(1 - cos 2ks) (24) 
in the region of stability of alterating ground state: 71 < 0, —71 + 472 > 0. 

The case of phase-modulated ground state is more complicated due to the 
lack of invariance with respect to the time reversal. The result in this case 
reads 

^3 — 7i(cos^o — COS (^0 + ks)) + 72 (cos 2^0 — cos2(^o + ks)) (25) 
and is valid in the third region of stability of ground states: 

72 <0, |7i|<4|72| (26) 
The quantity 6*0 in (25) equals the value of the variable 9 determined by (15). 
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7 Dispersion curves 



The calculation of the first and second derivatives of functions ei , £2 and 63 

describing excitation energy allows us to discuss shapes of dispersion curves 
in several sectors of the plane (71, 72) and to obtain the expressions for band- 
width W and effective mass rric- It is sufficient to plot the curves on the 
interval [0,7r]. The Fig. 1 contains examples of all 6 possible shapes of dis- 
persion curves. 

The first and second derivatives of excitation energy are expressed by the 
same formulae for both the east and nord-east sectors (the region of stability 
of the uniform state) : 

e'l = s(7i + 472 cos ks) sin ks, e'( — 5^71 cos ks + 472(2 cos^ ks — 1). (27) 

The effective mass, expressed by e'({k = 0), is then for both E and NE sectors 
given by the same formula: 

rnc = ^(7i + 472)-'. (28) 

However, the expressions for bandwidth are different. In the E sector e^>0, 
the dispersion curve is monotonic, and the bandwidth W is expressed by the 
difference of values of ei at the ends of interval [0, tt]: 

W = 2-fi. (29) 

In the NE sector the dispersion curve reaches maximum at the point ks — 9, 

where 6 is determined by (15), {9 < tt/2 for NW sector and 9 > tt/2 for NE 
one). The expression for bandwidth calculated as 61(6*) — ei(0) has the form 

W^ = 7i + 272 + ^. (30) 
872 

In the region of stabihty of the alternate state, i.e. in the nord-western and 
western sectors, the formulae for the first and second derivatives of excitation 
energy 62 have the form (27) with 71 replaced by —71. Hence the substitution 
of I71I instead of 71 into (29) gives bandwidth for western sector; the same 
substitution into (30) leads to the bandwidth for nord-western sector. The 
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expression for effective mass for both NW and W sectors may be obtained 
from (28) by tlie same substitution again and lias the form 

me=^(|7i|+472)-^ (31) 

For the region of stabihty of the phase-modulated foundamental states, 
i.e. for both southern sectors, we have 

= s[7i + 472 cos(6' + ks)] sm{e + ks) (32) 

and 

e'^ = 52^71 cos(^ + ks) + 472(2 cos^(^ + ks) - 1) (33) 

with 9 expressed by (15). In the interval [0,7r] €3 = if ^ + A;s = or 
6 + ks = TT, or else cos(^ + ks) = cos{6). That implies that €3 has minimum 
at ks = e {9 > tt/2 for SW sector and 9 <tt/2 for SE one). 

The expression for bandwidth in SW sector is equal to the difference 
^si^) ~ ^3(0) and in SE sector 63(71) — €s{9). The common formula for the 
bandwidth has the form 

W^=|7i|+2|72| + ^- (34) 

o|72| 



while the effective mass in both southern sectors is given by 

_ 4I72I 

S2 IQ^2 _ ^2 ■ 



(35) 



This completes the discussion of bandwidths and effective masses for all the 
types of ground states. 
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Fig. 1. Classification of dispersion curves, grade K=2 



